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The Abraham-Minkowski momentum controversy is the outwardly visible symptom of an incon- 
sistency in the use of the energy-momentum tensor in the case of a plane quasimonochromatic field 
in a simple linear dielectric. We show that the Gordon form of the electromagnetic momentum 
is conserved in a thermodynamically closed system. We regard conservation of the components 
of the four-momentum in a thermodynamically closed system as a fundamental property of the 
energy-momentum tensor. Then the first row and column of the energy-momentum tensor is pop- 
ulated by the electromagnetic energy density and the Gordon momentum density. We derive new 
electromagnetic continuity equations for the electromagnetic energy and momentum that are based 
on the Gordon momentum density. These continuity equations can be represented in the energy- 
momentum tensor using a material four-divergence operator in which temporal differentiation is 
performed with respect to ct/n. 

I. INTRODUCTION 

The energy-momentum tensor is a concise way to represent the conservation properties of a flow field. Simple 
in concept, the form of the energy-momentum tensor for the field in a dielectric has been at the center of the 
century-long Abraham-Minkowski momentum controversy [TJI2]- The tensor that was proposed by Minkowski [3] 
in 1908 is not symmetric and consequently does not conserve angular momentum [3]. Abraham [5] subsequently 
proposed a symmetric tensor at the expense of a phenomenological force. The original point of contention of 
the Abraham-Minkowski momentum controversy was whether the Abraham momentum density 

gA = (ExH)/c (1.1) 

or the Minkowski momentum density 

gM = (DxB)/c (1.2) 

provides the correct description of the momentum of light in a linear medium. The absence of an experimental 
decision allowed the debate to persist until the 1960s, when resolution of the Abraham-Minkowski dilemma was 
provided by Penfield and Haus [5] , based on earlier work by M0ller [7] , who showed that the issue is undecidable 
because neither the spatial integral of the Abraham momentum density nor the spatial integral of the Minkowski 
momentum density is the total momentum of the closed system. Likewise, the total energy-momentum tensor of 
the closed system is not the Abraham energy-momentum tensor and it is not the Minkowski energy-momentum 
tensor. 

What, then, is the total momentum and the total energy-momentum tensor? Here, we adopt the definition 
of the total momentum as the momentum quantity that is conserved in a thermodynamically closed system. 
This approach has the advantage of obtaining the total momentum directly from global conservation principles 
and neatly avoids any issues regarding the ill-defined roles of the Abraham and Minkowski momentums. We 
construct a thermodynamically closed system consisting of a homogeneous dielectric block illuminated by a 
quasimonochromatic field at normal incidence through an antircflcction-coating. In this closed system, the total 
energy and total momentum are conserved quantities and are well-defined by virtue of being conserved. We 
find that the total energy is the spatial integral of the electromagnetic energy density 

Pe = i(n 2 E 2 +B 2 ) (1.3) 

and that the total momentum is the spatial integral of the Gordon [8 momentum density 



g G = (nE x B) /c. 



(1.4) 
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The four-momentum density is therefore g a = (p e /c, gc)- We then show that the elements of the first row 
of the energy-momentum tensor can be chosen to be energy and momentum densities that satisfy the re- 
quirement for conservation of the components of the four-momentum. By doing so, the four-divergence of the 
energy-momentum tensor produces an energy continuity equation that is incompatible with Poynting's theorem. 
Retaining the conservation properties of the total energy-momentum tensor, we derive new electromagnetic 
continuity equations. We define a material four-divergence operator in which the temporal differentiation is 
performed with respect to ct/n and obtain a traceless symmetric energy-momentum tensor. The first-row and 
first-column components of the tensor are the densities of conserved energy and momentum quantities and the 
material four-divergence of the tensor generates the new electromagnetic continuity equations in terms of the 
electromagnetic energy density and the Gordon total momentum density. 



II. THE TOTAL ENERGY AND TOTAL MOMENTUM 



The first step toward obtaining the total energy and total momentum is to define a thermodynamically closed 
free-body system. We consider a stationary dielectric block illuminated at normal incidence from vacuum by a 
plane quasimonochromatic field. The dielectric has a gradient-index antireflection coating that allows radiation 
pressure on the vacuum/dielectric interface to be neglected. Then the refractive index n for a linear, isotropic, 
transparent dispersion- negligible dielectric is real and time- independent. The dielectric is homogeneous and 
occupies a finite region of 3-dimensional space so spatial variation of the refractive index, n = n(r), is limited to 
the transition region of the coating. The quasimonochromatic electromagnetic field is represented by the vector 
potential 

A(r,t) = - (A(r,t)e- l{uldt - kd - r) + k'faQe't"**-**")) (2.1) 
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where A is a slowly varying function of r and t, LOd is the center frequency of the field, = (nwd/c)ek, 
and is a unit vector in the direction of propagation. Figure 1 shows the envelope of the incident field 
Ai(z) — (A(z,£o) • A*(z,io)) 1 ^ 2 about to enter the dielectric through a gradient-index antireflection coating. 
Figure 2 shows a time-domain numerical solution of the wave equation at a later time when the refracted field 
A r (z) — (A(z, ti) ■ A* (z, ti)) 1 / 2 is entirely inside the dielectric. As shown in Fig. 2, the amplitude of the refracted 
field is A r = Ail y/n and the spatial extent of the refracted pulse is w r — Wi/n in terms of the width Wi of the 
incident pulse. It can be shown that this linear result is quite general in terms of the refractive index, as well 
as the width and amplitude of the field. 

For a stationary dielectric, the electromagnetic energy is 

U = [ ]- (n 2 E 2 + B 2 ) dv, (2.2) 

Jit 2 

where the integration extends over all space a. For the example shown in Figs. 1 and 2, there is no significant 
field outside the rectangular pulse. Then, the incident energy at t = to, 

U(t ) = ^A 2 , (2.3) 



2c 2 



is equal to the refracted energy at t — t\, 



U( tl ) = n ^A 2 r = ( 4 ) =R/(1 , (lUi 



>n 



The field has been averaged on a scale that is long compared to an optical period, but short compared to t\ — to, 
accounting for a factor of one-half. It can be demonstrated, in a similar manner, that the electromagnetic energy 
at any time t > to is equal to the incident energy. The temporal invariance of the electromagnetic energy makes 



U, Eq. (2.2), the conserved energy of the closed system. 
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Conservation of the electromagnetic energy is all that is needed to show that the Gordon momentum 

G G = / ^^dv, (2.5) 



obtained by spatially integrating the Gordon momentum density [8], Eq. (1.4 1, is the total momentum of our 



closed system. For monochromatic radiation, where B = nE, we can write the momentum as 

Gq = f —e k dv = (\ (n 2 E 2 + B 2 ) ^dv = -e k (2.6) 

J a C i/ cr c c 

in the direction of the propagation unit vector ek- If the total energy U is temporally invariant, and therefore 
conserved, then so is Gg = (U/c)ek. Alternately, we can show that the momentum balance 

G (*x) = ^A>e k = (*)\ = GoW (2.7) 



'ii 



is temporally invariant [9]. The momentum formula, Eq. (2.5), was originally derived in 1973 by Gordon [8]. 



Although there are some issues with Gordon's derivation, temporal invariance is decisive. Hence, the total 



momentum in our closed system is given by Gordon's formula, Eq. (2.5), and the total momentum density is 



the Gordon momentum density gc, Eq. (1.4) 9J. 



III. THE ENERGY-MOMENTUM TENSOR 

The energy and momentum continuity equations of a thermodynamically closed system can be combined to 
form a tensor differential equation and the energy-momentum tensor is the central clement of this formalism. 
While the energy-momentum tensor formalism is straightforward, it has not been successful in application to 
classical continuum electrodynamics. Now that we have identified the total momentum by global conservation 
principles [9] , we can use the apparatus of energy-momentum tensor theory to derive the energy-momentum 
tensor. However, it is not that simple. The Abraham-Minkowski momentum controversy arises from an incom- 
patibility between two of the properties of the energy-momentum tensor. 

The energy-momentum tensor has a number of important properties. Here, we employ the summation 
convention where identical indices on the same side of the equation are summed over, Greek indices belong to 
(0, 1, 2, 3), and Roman indices run from 1 to 3. Then, the four main properties of the energy-momentum tensor 
T a $ are: i) The four-divergence operator d a = (c~~ 1 d t ,d x ,d y ,d z ) applied to the rows of the tensor generates 
continuity equations 

d a T af} = (3.1) 

for the electromagnetic energy and the components of the momentum, ii) Conservation of angular momentum 
requires diagonal symmetry 

iii) The array has a vanishing trace 

T" = (3.3) 
corresponding to massless particles [H 110] . iv) The components of the four-momentum 

U= [ dvT m (3.4a) 



G i = - I dvT w (3.4b) 



c 
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are conserved [3] in an unimpeded flow. Conditions iv) and ii) dictate the elements of the first row and first 
column of the energy-momentum tensor, such that 



Pe cg Gl cg G2 cg G3 

ego i W n W 12 W 13 

cgc 2 W21 W 22 W 23 

cgoz W31 W32 W 33 . 



(3.5) 



where the elements of W are yet to be specified. Applying condition i) to the first row of the array in Eq. (3.5 1, 
we have 



ldpe 
c dt 

which is incompatible with the Poynting theorem 

dp e = 
dt 



V • (nE x B) 



V • c(E x H) . 



Alternatively, we can use Poynting's theorem and the momentum continuity equation 
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dt 



(D x B)/c + V • W = 



(3.6) 



(3.7) 



(3.8) 



to populate the tensor 



± M — 



Pe cg Al cg A2 cg A3 

CffMi Wn W12 W 13 

cg M2 W21 W 22 W23 

cgua w si w 32 W33. 



(3.9) 



where 



-EiDj - BiBj + - (E • D + B • B) 5 iS 



(3.10) 



is the Maxwell stress tensor. The resulting Minkowski energy-momentum tensor violates condition iv), in 
addition to condition ii). Because we are in a regime of negligible dispersion, we can rewrite the momentum 
continuity equation, Eq. (|3.8[), as 



d_ 
dt 



(E x H)/c + V ■ W = (1 



8 
] 0t 



(E x B) jc 



(3.11) 



and obtain the Abraham energy-momentum tensor, 



T 



a/3 



Pe eg A 1 cg A2 cg A3 

cg Al Wn W12 W13 

cg A2 W21 W 22 W23 

P9 A3 W31 W32 W33. 



Our condition i) becomes 



a/3 _ fA 



J a 



(3-12) 



(3.13) 



However, the Abraham force f£ = (0, c(l — n 2 )dtgA) is a source or sink of momentum such that momentum is 
not conserved in a closed system, also in violation of condition iv). 

We have demonstrated an inconsistency in the energy-momentum tensor formulation of classical continuum 
electrodynamics. Specifically we have demonstrated that all of the properties of the energy-momentum tensor 
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cannot be simultaneously satisfied for electrodynamic fields in a dielectric, even if the system is thermody- 
namically closed. The elements of the first row of the energy-momentum tensor can be chosen such that the 
four-divergence of the first row is the Poynting theorem or they can be chosen to satisfy conservation of the 
components of the four-momentum, but not both. 

In order to fully resolve the Abraham-Minkowski controversy, we must make the properties of the energy- 
momentum tensor self-consistent when applied to electromagnetic radiation in matter. Conservation of the 
momentum four-vector is paramount for the total energy-momentum tensor of a thermodynamically closed 
system. Likewise, conservation of angular momentum is required and we affirm our tensor properties iv) and 
ii) . Then the first row and column of the tensor are populated by the densities of the electromagnetic energy 
density and the Gordon total momentum density as shown in Eq. (3.5 1. These densities, integrated over space, 
correspond to the conserved total quantities of energy and momentum as shown in Eq. (3.4). Property hi) is 
not in question. 

The problem that remains is a situation in which the four-divergence of the total energy-momentum tensor, 
property i), is not consistent with the electromagnetic continuity equations, and Poynting's theorem in particular. 
We can multiply Poynting's theorem by n(r) and use a vector identity to commute the refractive index with 
the divergence operator to obtain a continuity equation 



n dp e 
c dt 



V • (nE x B) 



Vn 



(nE x B) 



The second-order energy continuity equation, Eq. (3.14), can be written as two first-order equations 



n 8{nE) 
c dt 



(V x B) 



(3.14) 



(3.15a) 



n dH _ N Vn , _. 

- — = - V x nE + x (nE). 

c at n 



(3.15b) 



Then, for a homogeneous dielectric, we can drop the term Vn x E and combine Eqs. (3.15 1 to obtain the energy 
and momentum continuity equations 



n d 
~cdt 



(n 2 E 2 



V • (nE x B) = 



(3.16a) 



77 (J 

-— (nE x B) + V ■ W = 0. 

c at 



The elements of W are the elements of the stress tensor 



(3.16b) 



Wij = -nE.nEj - B t B 2 + - (nE • nE + B • B) 5 ir 



We define a material four-divergence operator 



and replace property i) with 



afi _ 







(3.17) 



(3.18) 



(3.19) 



that generates Eqs. (3.16). Then the tensor, Eq. (3.5) is a traceless, symmetric energy momentum tensor whose 
material four-divergence generates the new electromagnetic energy and momentum continuity equations. It 
should be noted that the new energy continuity equation, Eq. (3.16a), is mathematically equivalent to Poynting's 



theorem, Eq. (3.7) because the two are related by a vector identity. Likewise the new momentum continuity 
equation, Eq. (3.16b), is mathematically equivalent to the momentum continuity equation, Eq. (3.8) if the 



radiation is sufficiently monochromatic and sufficiently far from any material resonances that dispersion can be 
neglected. 
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IV. SUMMARY 



A transparent linear dielectric block in free space illuminated by a quasimonochromatic field can be configured 
as an isolated system. However, the partial reflection of the field at the surface causes a change in the momentum 
of the field and the resulting radiation pressure accelerates the block. In principle, we could write a complete set 
of equations of motion at the microscopic level, but the effects of radiation pressure are not treated in a complete 
way at the level of the macroscopic Maxwell equations. In order to derive a theory of minimum complexity, 
we assumed a simple dielectric defined to be linear, isotropic, homogeneous, transparent, and dispersionless 
with negligible clcctrostrictive and magnetostrictive effects. The block is stationary in free space and radiation 
pressure on the antireflection coating is negligible. For this system, the elements of the energy-momentum 
tensor 



rpa/3 



Pe cg Gl cg G2 cg G3 

cgci W n W12 W 13 

c#g 2 W21 W22 W 23 

cga 3 w si W 32 W 33 . 



(4.1) 



are the electromagnetic energy dens ity p„ = (n E + B )/2, the Gordon momentum density go = (nE x B)/c 



and the Maxwell stress tensor, Eq. (3.101, in the form W{j = —nEinEj — B, L Bj + | (n 2 E • E 



properties of the energy-momentum tensor are: i) The material four-divergence operator 

- n d d d d 
c d t ' dx dy d z 

applied to the rows of the tensor generates continuity equations 

d B T afj = 



B B) Sij. The 



(4.2) 



(4.3) 



for the electromagnetic energy and momentum. The energy continuity equation and the momentum continuity 
equation are mathematically equivalent to Poynting's theorem and the momentum continuity equation for a 
linear, isotropic, homogeneous, transparent, dispersionless dielectric, ii) Conservation of angular momentum 
requires diagonal symmetry 



iii) The array has a vanishing trace 

corresponding to massless particles [H 1 10] . iv) The components of the four-momentum 

U = [ dvT 00 



(4.4) 



(4.5) 



(4.6a) 



c 

are conserved 0] in an unimpeded flow. 



G l = - / dvT 0t (4.6b) 
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Fig. 1. A quasimonochromatic field incidenton a linear 
medium with a gradient-index an ti reflection coating. The 
shaded region in the profile of the index of refraction. 
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Fig. 2. The quasimonochromatic field ofFig.l after it has 

propagated into the linear medium. 



